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A quantum mechanical close coupling theory of spin relaxation in the stretched F — 2, Mf = 2, 
hyperfine level of sodium is presented. We calculate the dipolar relaxation rate of magnetically 
trapped cold sodium atoms in the magnetic field range < B < 4 Tesla. The influence of shape 
resonances and the anisotropy of the dipolar interaction on the collision dynamics are explored. We 
examine the sensitivity of the calculated cross sections on the choice of asymptotic atomic state basis. 
We calculate and compare elastic scattering with dipolar relaxation cross sections for temperatures 
ranging from the ultra-cold to 2 K. We find that the value for the ratio of elastic to inelastic cross 
sections favor application of proposed buffer gas cooling and loading schemes. 
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I. INTRODUCTION 

Advances in the cooling and trapping of atoms have greatly facilitated the exploration of quantum degenerate 
matter nj. The realization of Bose-Einstein condensation (BEC) |_l y. |j. _f in atomic vapors validates the standard 
theory 6] for non-interacting and weakly interacting systems, but experiments 0, _1 13 demonstrate that atomic 
interactions, though weak in an ensemble of atoms in the gas phase, lead to interesting phenomena that are not 
present in the ideal gas system 0, 0, m> 0] . 

In a dilute, cold, gas, atoms interact primarily via long range dispersion and exchange forces. However, inelastic 
processes are driven by spin exchange 0, __J and dipolar interactions . The latter process does not conserve total, 
atom pair, spin angular momentum and it is a primary mechanism by which atoms, having hyperfine structure, can 
suffer an inelastic transition. Dipolar relaxation determines the lifetime of the hydrogen atom BEC [l9l | , contributes 
to heating and influences the o per ation of atomic clocks |2d| . Rates for dipolar relaxation have been measured in 
7 Li [2l|, Cs [2^, 85 Rb 01 j H |19| and Cr [2j|. The rates are generally small, typically having values that range 
10 -14 — 10 -16 cm 3 s _1 , but in the cases of Cs and Cr anomalously large values have been reported [22I I23L l24j . 
Calculations for dipolar relaxation rates, in the zero temperature limit, of several species have also been reported 

m mm mm mm _____ 

New magnetic trapping and buffer gas cooling schemes |2a, __} create the opportunity to study a host of atomic 
and molecular species that are not amenable to laser cooling technology. In a typical loading scheme, species with 
large magnetic moments are trapped by external fields at relatively high temperatures, on the order of 1 K, before 
they are cooled into the sub-Kelvin regime. In order to model this process one needs a detailed understanding of 
the collision processes that can lead to trap loss and heating. To that end, we present a comprehensive quantum 
mechanical theory of dipolar relaxation in alkali atoms. The theory is suited for application in gases at temperatures 
where many partial waves in the collision wave function contribute, and is applicable at arbitrary external magnetic 
field intensity. We apply the theory to calculate the dipolar relaxation rate of the stretched hyperfine level in the 
23 Na(3s) atom and, in this paper, present results for temperatures that range from the ultra-cold to several Kelvin 
and magnetic field intensities in the range < B < 4 Tesla. The results of our calculation are compared to previous 
theoretical predictions |29l l3l). We present, the first fully quantum mechanical calculation for dipolar relaxation of 
sodium atoms in a magnetic field at higher temperatures where many partial waves contribute. We give a detailed 
description of the collision theory and explore the consequences of anisotropy on the collision dynamics. We point 
out the importance of shape resonances and their influence on the value of the inelastic rate. We identify a feature 
in the cross section that corresponds to the presence of an above threshold resonance, or a virtual state, in the I = 2 
partial wave of the scattering amplitude. 

In section I we provide an introduction to the theoretical formalism that is applied in the calculations. A detailed 
discussion of the close coupling equations, asymptotic boundary conditions and symmetry requirements is given in 
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TABLE I: Quantum numbers associated with the various basis representations. Mf is the total spin angular momentum along 
the quantization axis and Ay is the energy defect between the F = 2 and F — 1 hyperfine levels in Sodium. 
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sections II, and III. In section IV we present the results of our calculations, and provide a detailed analysis of these 
results. Unless it is otherwise stated, atomic units are used throughout the discussion. 

II. CHANNEL BASIS 

We consider two sodium atoms in their F = 2, Mp = 2 hyperfine level, the maximal stretched state. In Table I we 
itemize two-atom hyperfine levels with the notation \FMpf a f b >, where F is the total angular momentum of the two 
atoms, Mp the azimuthal projection of that angular momentum, and f a , f b are the total angular momenta for atoms 
a and b respectively. The states listed can mix, through dipolar and spin exchange interactions, with the maximal 
stretched state during a collision. 

Dipolar interaction selection rules(discussed in the sections below), allow a change in the azimuthal quantum number 
Am f = 2, 1 and thus the states itemized in Table I must be included in the close coupling expansion. States within 
a given Mp manifold can undergo spin-exchange transitions. In Table I we also list the basis Ifafnafbin-b > in which 
the individual atom azimuthal angular momenta are good quantum numbers. The basis \SMsIMi > diagonalizcs 
the asymptotic Hamiltonian if the hyperfine interaction can be neglected, i.e. at large magnetic field strengths. Here, 
S, Ms are the total two-atom spin angular momentum quantum numbers, and /, Mj the nuclear angular momentum 
quantum numbers. Allowed values for 23 Na2 are 5 = 1,0 and / = 3, 2, 1, 0. In the close coupling expansion involving 
molecular channel states we keep the notation that is appropriate for the asymptotic region to itemize the states 
in the expansion. For example, we define molecular channel basis \SMsIMi >= | 3 S„ > ®\IMj > for S = 1, 
and \SM s IMi >= ^Ep > ®|/Mj > for S = 0, where | 3 S„ >,| 1 S g > are Born-Oppenheimer (BO) eigenstates for 
the ground Na2 system. The molecular channel basis merge to the correct asymptotic basis at large inter-nuclear 
separation. The states \FMpf a f b > are then defined by the linear combination of the BO channel states given above, 

\FM F f a f b >= \SMsIMj >< MjIMsS\FM F f a f b > (1) 

where the coefficients < MiIMsS\FMpf a f b > are standard recoupling coefficients appropriate for the asymptotic 
basis. In this way we insure the molecular close coupling expansion accounts for the asymptotic hyperfine interaction 
within each atom. In the case of a homonuclear system, the basis \FMpf a f b > is not an eigenstate of the electron 
inversion operator, but we can define the states \FMp(f a f b ) >= -^{\FMpf a f b > ±\FMpf b f a >, that are eigenstates. 

In a non-zero magnetic field, the asymptotic Hamiltonian is not diagonal in the representation defined by the 
basis vectors introduced above. In addition to asymptotic hyperfine interactions, each atom experiences the Zeeman 
interaction. In that case, a linear combination of the basis states defined above must be found so that the asymptotic 
Hamiltonian is diagonal in the new representation. We express these states using the notation \Mppei >, where 
Mp is the total angular momentum along the direction fixed by the magnetic field, p is an inversion parity quantum 
number, and is the asymptotic energy level eigenvalue. We discuss the construction of that basis in the sections 
below. 

In Figs, la, lb we illustrate the energy spectrum of the asymptotic Hamiltonian for the Na2 system, as a function 
of magnetic field strength. 



3 





B(Tesla) 



FIG. 1: (a) Energy levels for a pair of sodium atoms in a magnetic field. The levels correspond to the states itemized in Table 
II. In the figure they are grouped in order of decreasing energy, (b) In the large B field limit the arrows represent the total 
electronic, azimuthal, spin angular momentum of the corresponding levels shown. 



III. CLOSE COUPLING EXPANSION 



In the Pauli approximation, the magnetic Breit interaction between the two valence electrons is given by [33j 
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where a is the fine structure constant, Si the spin of electron i and the displacement vector for the two electrons. 



In order to include the magnetic interactions in the scattering equations we substitute r; 



R, where R is the 



inter-nuclear vector of the two atoms. This approximation is valid at large inter-nuclear separations and we replace 



4 



expression by the model interaction 

,(Si-R)(S 2 -R) 



Hdi— polar — 



Si-Sa-3- 
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where we have ignored the Fermi-contact term. We did not include the electron spin-nuclear spin and the nuclear 
spin-nuclear spin interactions since they contribute to the interaction energy an amount that is at least three orders 
of magnitude smaller than the interaction energy obtained from Eq. ffifr. 

Using standard Racah- algebra techniques, we re-express Hdi- po iar in terms of irreducible tensor operators, thus 


H dl - polar = v(R) ^(-i)^y( 2 ) ? (^)5^ (4) 

9 

where Y^ 2 ' are the components of the spherical harmonic of rank 2, is a second rank tensor in the product space 

spanned by the states \S1m1S2m2 >, and v(R) = expressed in atomic units. 

In the special of case of a null external magnetic field, both \FMpf a f a > and \ f a m a fbmb > basis vectors, itemized 
in Table I, diagonalize the total Hamiltonian in the asymptotic region. We use the former set to express the system 
wavefunction by the close coupling expansion, 

*(R, r) = F FMF f afb (R)\FM F f a f b > (5) 

where the sum is over all quantum numbers itemized in Table I, and from which we obtain the coupled equations, 

- ^V 2 F,(R) + J2 ^m( R ) + E Kf^(R) + E K7( R )^( R ) = EFi(H). (6) 



In deriving Eq.© we have ignored non-adiabatic effects[l3,[T6(, since they are expected to provide a small correction 
in the sodium-sodium system. The subscript on the scattering amplitude denotes the channel quantum numbers, 
i = {FMpfafb)i, E is the total energy, and V^(R), V^(R), Vy*(R) are multi-channel potentials that correspond 
to the electrostatic, nuclear hyperfine, and di-polar magnetic interaction Hamiltonian respectively. The explicit 
electrostatic and hyperfine terms are 

V(*) = J2(-l) fa+fb+f '« +fl Sf,F> S MfMf , [F, S, I][faJa'J b , X 
SI 

'1/2 3/2 f a ) f 1/2 3/2 /„, 

1/2 3/2 f b } { 1/2 3/2 f v ) e s (R) (7) 
S I F 5 / F 



V i3 Hf = [/„(/« + 1) + h{f b + 1)] (8) 

where A F is the Fermi hyperfine constant for the ground state of sodium, /1 is the nuclear reduced mass flH llfij . 
and es{R) are the Na 2 ground state Born-Oppenheimer potentials, for the triplet 5=1, and singlet S — states 
respectively. In deriving Eq. J7J we used the relation 

<IM I SM s \FM F f a f b > = [F,f a J b ,S,I}^ 2 (-l) F+M " x 

. X f 1/2 3/2 f a ) 

We now derive an expression for the di-polar interaction. Using Eq. W e have 

U™(R) = v(R) 53("1)% (2) < fa-fvM F ,F'\S [ *l\FM F f a f b > (10) 
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We express amplitude JHJ by a partial wave expansion in spherical harmonics, 



lm 



R 



(14) 



(15) 



Though Vfj, V t j are isotropic in the nuclear orientation Vy 1 , according to Eq. i|14(l . is not and the partial wave 
expansion does not lead to radial equations that are diagonal in the nuclear angular momentum / and m. Inserting 
(HSJ into Eq. © we obtain, 



E E Uij{lm,l'm' ',R)F j (l'm',R) = EF l (lm,R) 

j V ,m f 
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According to Eq. i|14[l Q = Mf — Mf> = AMj?, and from Eq. (|17|l q = rn — m! = Am. Therefore we obtain the 
selection rule Am = AMp and from the selection rules for the 3j symbols in Eq. l|ltj|) we require I — V = Al = 0, 2 , 
and 1 = 1' ^0. 



IV. SCATTERING FORMALISM 



It is useful to re-express the coupled radial equations in the form 

+ T(R)E(R) + V hf F(R) + v m (i?)F(i?) 



i . d 2 1(1 + 1), 



2y, K dR 2 R 2 

EF{R). (18) 

In the notation introduced above F(R) is a square matrix whose columns contain the independent solution vectors to 
the coupled equations (|16fl . The row and column indices for matrix F(i?) itemize both the internal and orbital angular 
momentum quantum numbers. A given value of index i, identifies the set i = {(FMpfalbjihrni} where li, uij are 
the total and azimuthal quantum numbers for channel i. At a given collision energy, we are allowed to truncate the 
partial wave expansion (|15(l at some maximum value l m ax and matrix F(i?) is a finite n-dimensional square matrix 
where n = rii \\_{^ (21 + 1), and rii is the dimension of the internal Hilbcrt space. 

Matrices V e (i?), V , correspond to the electrostatic and hyperfme Hamiltonian respectively, and they are diagonal 
with respect to the angular momentum quantum numbers. The matrix 1 is diagonal and contains the channel angular 
momenta along the diagonal. However, V m (i?), whose components are V™(i?) = Uijilim^ Ijirij; R), is not diagonal. 

In the limit R — > oo we require that 

F l? (i?) -» -^=\ S. ir sin(hR-h^)+K llC os(hR-h^) \ (19) 



where K- arc the elements of the K matrix. 
— l j 

We introduce the amplitude G(R) = F_(R) C where C is a constant matrix chosen so that in the limit R — > oo 

Gy(JJ) - ^|*yrap(-i(^fl-l i |))-S w ejcp(i(fc i iJ-/ i |))| (20) 
where the radial S-matrix S is, 

S = Q-tKJG + iK) -1 (21) 

We construct a reduced multichannel amplitude, Gryi(R), where the notation [ij] implies that the indices denote the 
quantum numbers of the internal states only. We define 

%]( R ) = E E *W9fl*WPifr) ,1/2 R ( 22 ) 
and find that in the asymptotic limit R —> oo 

G fo] (R) - exp( l K. i ■ R) + f m (00; 0,0,)^^ (23) 

where we have used Eq. (|20|l . and /r^-i (00; 0i0j) is the scattering amplitude for the system to undergo a transition 
from an initial internal state j into an internal state i and into solid angle d6(sin6)d<j) following an initial approach 
along the incident wave vector Ki with polar angles 0j0j. Comparing expression ill'L'l) and (|23J) we find that 

f m (00; 4 2 ) = J2 E (^)^i^(^ ~ (24) 

Though Gryi(R) has the desired asymptotic behavior for scattering solutions it does not posses the symmetry 
required by the Pauli principle. Because the sodium nuclei are identical fermions, the total wavefunction must be odd 
under their interchange. Let P^ 2 be the nuclear permutation operator, thenjl8l| 

P»\FMf a ft >= (-l) F+f ° +h+1 \FMf b f a > . (25) 
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We introduce a shorthand notation for the channel indices that label the matrix G(R); if i = FMf a fb then i = 
FMftfa- In this notation the above relation is written P^\i >= (— l) i?+ ^ a+ ^ 6+1 |i >. If the system is initially 
prepared in state j, and is given by I* >= J2i > i then we require that R) = (— l^+^+^G^-i (R). 

Using this notation we replace (1221) with, 

liTTli IjTTlj j 

| G^R) + { _ 1)F+h+fa+fb Gy{R) 1 

which has the desired symmetry. 
In the asymptotic limit we get, 

G fo] (R) - 6 M exp( i K l -H) + (-l) F +^+^5 m cxp(-iK l -Il) + 

fmQMik) + {-l) F+U+fb f m {-e + + tt; w] exp ^ (27) 

where we have used fcj = fcj. The cross section for a system in an internal state |j >= \F'M' f a i fv > to undergo a 
transition into state |z >= \FMf a fb > is 

, . Vi 1 1 

F+f a +f b t J_o , _ J. , _. /W ,\|2 



/ dihlftolW-Afa) + (-1) (-6* + 7T + tt; 1 (28) 

where we integrate over all scattering angles and average over all directions of the incident wave. Vj is the velocity 
in the incoming channel and Vi the final channel velocity. We have included a factor of \ in order to insure that the 
incoming flux is normalized to unity. Using expression (|26|l we can re- write Eq. Q27[l. 

°(l "» *) = ^EE \ T M{Umi\ljrnj) + (-l) F+/ « +/i - + ^T fo . ] (/ i m i ;/ i m J -)| 2 

3 limi Ijmj 

T[ij]{kmi\ljmj) = S[ij] — S[ij] (29) 

We use Eq. (|29|l to calculate the total inelastic transition cross section in the case for zero, or small, magnetic field 
intensities. The initial state corresponds to the maximal extended state |F = 4M = 4/ a = 2/b = 2> and at low 
energies only incident s-waves contribute. According to the dipolar selection rules the exit channels are d-waves, and 
we obtain a simple expression for the total inelastic cross section 

a T = -^\^FM F f a f h {mi)\ 2 (30) 

FM F f a f b m 4 =-2 

where, f FMpfafb (m 4 ) = T [ij} (k = 2, m 4 ; lj = 0, rrij = 0) for i = FM F f a f b and j = F = 4 M F = 4, /„ = 2 f b = 2, and 
k = kj. In deriving Eq. we used the fact IV, = (-l) F +f°+f<>T [i:j] . 

For a large magnetic field, such that the Zeeman splitting is much greater than the hyperfine interaction, we 
construct close coupling equations by using the basis vectors \SMsIMj > in expansion Eq. ©. We obtain an 
equation analogous to Eq. I|18|l except that V_ h * is replaced by an expression that describes the Zeeman interaction 
with the external field. In addition, the electrostatic and dipolar interaction matrices are replaced by V_ (R) and 
V m (R) respectively. They are related by the unitary transformation V_ e (R) = UV e {R)lT\ V_ m (R) = UV m (R)lJ-\ 
where =< FMpfafbljmjlkmiSMsIMj >. 

Because the states \SMsIM/ > are eigenstates of the nuclear interchange operator P^, i.e. 

P&\SM 3 IMi >= (-l) s+I+1 \SM s IMj > (31) 

we obtain 

°(J -») = ^EE \T M (km t ;l jmj )(l + (-l)'+^)| a (32) 

3 liVTLi IjTTlj 
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TABLE II: Quantum numbers associated with the states \MfP^ > that diagonalize the asymptotic Hamiltonian, Eq. 1321 . 
We itemize only those states whose parity, under nuclear interchange, is odd. The parameter £ = 4 ^ F B where B is the magnetic 
field strength and Af the energy defect between the F = 2 and F = 1 hyperfine levels of Sodium. The last column itemizes 
the B ^ limit of the states expressed in the \FMp(fafb) > representation. 
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where the channel indices now itemize the states in Table 1 under the \SMsIMj > representation. 

If the magnetic interaction energy is of the same order as the hyperfine energy, neither the \FMpf a fi, > nor 
the \SMsIMj > representations constitute a valid asymptotic basis since off-diagonal elements persist at large 
inter-nuclear separations. Instead, we choose a linear combination of these states that diagonalize the asymptotic 
Hamiltonian, 

H hf + H z 

H z = 2^ B B • S (33) 

where H h * is the hyperfine interaction, B is the external magnetic field whose orientation defines our lab quantization 
axis, S is the total electronic spin for the atom-atom system and \ib is the Bohr magneton. We ignored the magnetic- 
nuclear term since it provides a considerable smaller contribution to the total magnetic interaction energy than that 
given by Eq. I|33|) . The diagonalization procedure can be carried out numerically, and in Table II we itemize those 
states which contribute to dipolar loss from the incident, extended state, channel. Good quantum numbers for these 
states include the total azimuthal quantum number Mp, the nuclear interchange parity, and the energy eigenvalues 
itemized in Table II. The extended state is odd under nuclear interchange i.e., P^\Mp = 4; p = 1 e\ >= (— l) p \Mp = 
4; p = 1 ei > where we have used the notation described in Table II. 

Only states of odd parity are allowed as exit channels and these states are listed in Table II. Invoking the procedure 
discussed above, we obtain for the total dipolar loss cross section 

m i= 2 

°*=E E F^ (m<)|2 (34) 

j mi=-2 

where the sum over index j denotes the channels itemized in Table II. 
The rate coefficient for dipolar relaxation is given by the expression 



where p, is the reduced mass of the 23 Na2 system and <jt{E) the total inelastic cross section expressed as a function 
of collision energy. 

V. RESULTS AND DISCUSSION 

In Fig. 2 we plot the total rate coefficient (solid line), in the T — > limit, as a function of magnetic field strength. 
In Fig. 2 we notice that the total relaxation rate is nearly constant for field strengths up to about 100 Gauss(G). In 
the range 100 G < B < 400 G the rates exhibit significant structure. For larger values of B the total rate diminishes 
in a monotonic manner. In figure 2 we also plot, shown by the dashed line, results obtained using the approximation 
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FIG. 2: Total dipolar relaxation rate (heavy solid line) as a function of external magnetic field strength. The dashed lines 
corresponds to the approximation where the hyperfine basis \F Mf fafb > is used for the asymptotic channel states. The dotted 
line corresponds to the case where the \SMsIMi > basis is used and the hyperfine interaction is ignored. 



where the hyperfine states \FMp f a fb > define the asymptotic basis states and the asymptotic off-diagonal terms, due 
to magnetic interactions, are neglected. For small B the approximation gives excellent agreement, for the total dipolar 
loss rate, with the results obtained using the appropriate \MfP£ > basis. However, for B > 60 G this approximation 
considerably overestimates the total rate. 

The dotted line in the figure presents the results of the calculation when we used the \SMsIMj > representation to 
define the asymptotic channel basis, and where we ignored the hyperfine Hamiltonian. For small B this approximation 
is poor because it neglects the dominant hyperfine interaction. In the B — > limit, rates obtained in this approximation 
vanish due to the nature of the anisotropic dipolar interaction. According to the selection rules discussed in the 
previous sections, a pair of atoms approaching as an s-wave must exit as a d-wave, and therefore, the exit channel 
must be exothermic with respect to the entrance channel for the collision to proceed in the zero temperature limit. 
At B — > a finite energy defect is generated by the hyperfine interaction and, if hyperfine effects are ignored, the 
dipolar rate vanishes in the T — > 0, limit. This effect is clearly evident in Fig. 2. As B — * oo the neglect of the 
hyperfine interaction is justified and, in that case, we expect that the rate obtained using the \SMgIMi > basis to 
be a good approximation. In Fig. 2 we note that for B > 1 T the dotted line merges with the solid line and illustrates 
the validity of that approximation at large field strengths. 

To understand the nature of the observed structures in the total rates we neglect the fine structure and study the 
collision dynamics in the \SMsIMi > basis and show the results in Fig. 3. In that figure, the dashed line denotes 
the partial rate into the state where both atoms flip their total electronic spin, whereas the solid line corresponds to 
the case where only one atom flips its spin. We first consider the kinematics of the latter case. 

In the \SMsIMj > basis the Zeeman energy splitting between the IS* = l;Ms = 1;J = 3; Mi > and the \S = 
1; Ms = 0;I = 3; Mj > exit channel is given by 



AE= \2^ B B\ 



2{i 



(36) 



where k is the wavenumber that corresponds to the kinetic energy of the system in the entrance channel and B is the 
absolute value of the magnetic field that is parallel to the laboratory z-axis. In the k — > limit the exit channel is 
a d-wave and, using the potential for the 3 S U state of 23 Na2 system tabulated by Samuelis et al. |3fi| . we find that 
the 1 = 2 centrifugal barrier has a height AE(l = 2) = 1.6585 x 10~ 8 a.u. We equate the Zeeman splitting with the 
barrier height and find that the critical magnetic field strength required so enough kinetic energy is available in the 
exit channel to overcome the barrier has the value B c = 39.0 G. 

According to Fig. 3, at this field strength the relaxation rate is rapidly increasing, as B increases, but it is in 
a region to the left of its maximum which occurs at B max = 46.0 G. Therefore, the pronounced structure seen in 
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FIG. 3: Resonance/threshold structures observed in the dipolar loss rates that were calculated ignoring hyperfine effects and 
in the \SMsIMi > basis. The dashed line corresponds to dipolar loss involving total electronic spin flips for a single atom, 
whereas the solid line corresponds to spin flips involving both atoms. 



this rate cannot be solely attributed to a threshold effect. Indeed we found a resonance in the I — 2 partial wave 
that is due to the existence of a virtual state at E v — 1.85 x 10~ 8 a. it. Using Eq. I|36|l to convert E v into a field 
strength, we find B* = 43.5 G, a value that is about midway between B c and B max . For transitions into the state 
\S = 1; Ms = — 1; I = 3; Mj >, whose rate coefficient is shown by the dashed line in Fig 4, B we evaluate B* = E v /4/iB 
and get B* = 21.8 G. This value is close to B max = 23.0 G seen in the figure. These observations strongly suggest 
that the structure, evident in the in Fig. 3, is a consequence of shape resonance phenomena [3(||33]. This conclusion 
is strengthened by studies, discussed below, of this collision process at higher temperatures. 

Our rate in the B — > limit is in harmony with that that reported in a previous study [31| , but several times larger 
from that predicted in an earlier study |2^. In Refs. p9l l3l|j the authors used the \ f a m a fbmb > basis to calculate the 
rates for B ^ 0. We have shown here, that this approximation overestimates the relaxation rates for B > 60 G. The 
largest uncertainty is probably associated with the choice for molecular potentials of the ground 23 Na system. We 
adopted the most recent, and accurate, potentials tabulated by Samuelis et al. |35| . 

In Fig. 5 we present the results of our calculation for both the elastic and inelastic, dipolar relaxation, cross sections 
in collisions of spin polarized sodium atoms. The collision energies considered range from ultra-cold to 2 Kelvin. At 
higher temperatures many partial waves contribute and the multi-channel, close coupling theory described in the 
previous sections is applied. We display results obtained for magnetic fields that range between 1 and 4 Tesla. In 
the previous paragraphs we justified the neglect of the hyperfine interaction, in the calculations for field strengths 
B > IT. The results shown in Fig. 5 are based on this approximation, but because of the anisotropy in the dipolar 
interaction, this calculation still involves a large number of channels since coupled partial wave angular momenta up 
to J ~ 20 are required for convergence at temperatures T > IK. We find that the elastic cross sections dominate and 
are largely insensitive to the value of the applied field. In the range of applied fields considered and in the fiK collision 
energy region, the dipolar cross sections are about 10 -3 smaller than the elastic cross sections. The ratio decreases 
at higher collision energies and at cryogenic temperatures the ratio of elastic to inelastic cross sections « 10~ 6 . At 
higher temperatures we find that the dipolar relaxation cross sections decrease as the applied field is increased. In 
Fig. 5, we note that the elastic cross sections tend to a constant value as the gas temperature approaches the /jK 
range whereas the relaxation cross sections increase, in conformity with the Wigner threshold laws. Both the elastic 
and inelastic cross sections display resonance features discussed in the previous paragraphs. Because of the favorable 
ratio of elastic to inelastic cross sections, spin polarized sodium is, potentially, a good candidate for buffer gas loading 
and evaporative cooling. 
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FIG. 4: Elastic and inelastic cross sections in collisions of Na(F = 2Mf — 2) + Na(_F = 2Mp — 2). The collision energy is 
expressed in units of Kelvin. 
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